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PHASE  ACCURACY  EXPERIMENTS  WITH  A 
DIRECT-SAMPLING  COHERENT  DETECTOR 


INTRODUCTION 

The  coherent  conversion  of  radar  and  communications  signals  into  digital  form  may  be  accom¬ 
plished  by  frequency  translating  the  signal  to  in-phase  and  quadrature  (/  and  Q)  baseband  channels, 
each  containing  an  analog-to-digital  (A/D)  converter,  or  by  sampling  and  digitizing  the  signal  directly 
in  the  radio  frequency  (RF)  or  intermediate  frequency  (IF)  sections  of  the  receiver.  Considerable 
interest  in  the  latter  approach  has  developed  in  recent  years  [1-5]  since  A/D  converters  have  become 
faster  and  less  hardware  is  needed;  direct-sampling  methods  have  been  shown  to  be  more  accurate, 
especially  if  signals  vary  over  a  wide  dynamic  range. 

This  report  describes  a  direct-sampling  coherent  detector  model.  The  model  is  capable  of  con¬ 
verting  samples  of  radar  or  communications  signals  having  intermediate  carrier  frequencies  up  to  30 
MHz  [1]  and  bandwidths  up  to  5  MHz  to  corresponding  digitized  samples  of  /  and  Q.  Tests 
described  here  demonstrate  phase  accuracies  better  than  0.3°  at  IF  frequencies  below  4  MHz. 

The  report  briefly  reviews  the  rules  of  direct  sampling  to  clarify  the  method  used  to  derive  /  and 
Q  values  from  these  samples.  Special  attention  is  devoted  to  the  general  problem  of  measuring  phase 
accuracy  since  the  estimated  accuracy  of  the  device  under  test  exceeds  that  of  readily  available  phase 
standards.  Test  instrumentation  and  results  are  presented,  followed  by  a  discussion  of  the  effect  of 
phase  noise  and  A/D  converter  aperture  uncertainty  on  the  experimental  data. 

SAMPLING  FUNDAMENTALS 

If  the  spectrum  R(f)  of  a  signal  r(t)  contains  no  energy  outside  a  band  from  -  W  to  +  W,  r(t) 
may  be  represented  without  loss  of  information  by  samples  r(tn),  taken  at  points  spaced  by  no 
more  than  1/2  W.  That  is,  r(t)  may  be  exactly  reconstructed  from  these  equally  spaced  sampled 
values,  r„  =  r(tn).  The  above  two  sentences  state,  in  terms  appropriate  for  this  report,  the  sampling 
theorem  given  by  Shannon  in  1949  [6] . 

A  similar  theorem  may  be  stated  for  bandpass  signals.  If  the  spectrum  R(f )  of  a  signal  r(t) 
contains  no  energy  outside  the  bands  specified  by  the  inequality  {M  -  1)  W  ~  |/|  <  MW  (M  is 
any  positive  integer),  then  r(t)  may  be  exactly  reconstructed  from  sampled  values  r„  taken  at  points 
spaced  by  r„  -  tn  ,  =  1/2  IT.  Note  that  the  sampling  rate  in  this  bandpass  case  must  equal  2  IT;  in 
the  lowpass  case,  2IT  is  simply  a  lower  bound  on  sampling  rate. 

It  is  worth  mentioning  in  passing  that  this  restriction  may  be  avoided  if  second-order  sampling* 
is  performed.  However,  this  involves  unequal  sample  spacing,  and  we  therefore  confine  our  attention 
to  the  case  where  energy  is  confined  to  the  bands  defined  bv  (M  -  1)IT  =  |/1  =  MW. 

Manuscript  approved  October  2X.  l^KX. 

*  Linden  s  paper  J7|  provides  a  lucid  review  of  sampling  theorems  Although  minimum  sampling  speed  is  often  called  the  Nwposf  utc. 
Shannon  fftj  is  usual))  credited  with  proving  that  a  band  limited  signal  may  he  exactly  reconstructed  from  samples 
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1  his  restriction  on  the  frequency  band  imposes  a  corresponding,  but  not  damaging,  restriction  on 
center  frequency .  /„  -  ( 2 A/  1)  U'/2.  Since  M  is  a  positive  integer,  the  allowed  values  of  /0  may 
,ve  written  IV  2.  31V' 2.  5VV  2  .  .  . .  Hence,  the  system  designer  may  select  an  IF  that  is  high  enough  to 
avoid  filtering  problems  but  low  enough  to  avoid  problems  with  A/D  converter  aperture  uncertainty. 
Comments  on  this  trade-off  are  included  in  the  final  section  of  the  report. 

DIGITAL  COHERENT  DETECTOR  (DCD)  DESIGN 

Figure  1  plots  a  segment  of  a  Gaussian  pulse-modulated  sine  wave  to  illustrate  the  locations  of 
direct-sample  points  when  M  -  3.  The  pulse  is  time-shifted  arbitrarily  so  that  tt  peaks  at  t  =  3:  its 
c  ‘trier  frequency  is  2.5  with  7t/4  phase  shift,  and  the  Nyquist  bandwddth  is  W  =  1. 


f  it;  I  -  Sine  wave  modulated  by  Gaussian  pulse.  Sample  points  are  labeled  with  corresponding  /.  Q  values. 


Figure  2  is  a  simplified  block  diagram  of  the  DCD.  The  actual  bandwidth  B  of  the  signal-plus- 
noise  is  determined  by  the  receiving  filter  and  must  be  less  than  W.  Digitized  sampled  values  (r„)  are 
processed  as  follows  to  provide  corresponding  /„  and  Qn  values. 

A  sample  point  is  selected  to  correspond  to  i  =  0.  Since  the  /„  —  ( —  I )'' r„  for  even  n  under 
these  assumptions  fl.2|,  /,,  =  rn,  I2  =  —  r2  ....  as  shown  in  Fig.  1.  When  n  is  odd. 
(J„  -  (  I)"'  u'~r„.  This  leads  to  the  Qn  values  noted  in  Fig.  1.  Hence,  for  these  /„  and  Q„ 
value  >  Q  | .  I2.  Q\  .  .  .).  the  only  computing  necessary  consists  of  switching  the  sampled  value  to 
the  correct  output  port. 

DCD  applications  often  demand  accurate  phase  where  the  values  Q{),  / 1 .  Q2,  />.  ...  must  be 
“■■ipuie-l.  I  h i s  computation  requires  a  digital  filtering  operation.  Another  investigator  1 2 )  has 
developed  a  solution  to  thi>>  problem  based  on  the  Hilbert  transform.  The  method  described  here  and 
m  Rets  I.  4.  and  5  uses  a  modification  of  the  sampling  theory  interpolation  function,  sometimes 
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r(n/2W) 


SAMPLING  PULSES 
Fig.  2  —  Basic  coherent  detector  approach 


called  a  frequency-window  function.  Both  methods  provide  exact  values  if  the  interpolation  is  based 
on  arbitrarily  long  data  sequences.  However,  practical  interpolation  must  be  computed  on  a  finite 
block  of  sample  data.  Thus  we  used  work  by  Helms  and  Thomas  |8]  based  on  what  they  term  “self- 
truncating"  modifications  of  the  conventional  sampling  function  (a  cardinal  function  in  the  low-pass 
case).  The  thrusts  of  their  work  were  the  improvement  of  interpolation  accuracy  and  the  reduction  of 
truncation  length,  hence  the  length  of  the  interpolation  filter  impulse  response. 

Figure  3  shows  the  means  used  to  compute  values  of  Qn  for  even  n  and  /„  for  odd  n.  We  con¬ 
tinue  to  deal  with  the  example  of  Fig.  1  (M  =  3).  By  using  finite  impulse  response  (FIR)  filter 
weights  (determined  by  the  above-mentioned  interpolation  function)  we  estimate  a  value  of 
r{tn  +  1/4  f'f).  If  /<;  »  W,  this  time  shift  of  1/4  /(J  is  equivalent  to  a  90°  phase  shift;  hence. 
r(0  +  1/4  f(\)  =  Qn.  Since  f\  is  a  fictitious  carrier  frequency,  it  can  be  made  arbitrarily  large;  this 
makes  (-l)"/2r(/„  +  1/4  f\)  arbitrarily  close  to  Qn  for  even  n.  Similarly. 

(-1)*''  "  >,/2rU„  +  1/4  /( J)  approaches  /„  for  odd  n  { 1 ,4J. 
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The  FIR  filter  shown  in  Fig.  3  computes  an  estimate  of  r(/  -+•  1/4  /, j)  based  on  samples  of 
r(n  It  happens  that  the  interpolation  function  is  odd  symmetric,  with  zeros  at  even-valued  displace¬ 
ments  from  its  center  (stages  2.  4.  6.  8.  and  10  in  Fig.  3).  Hence,  only  three  multiplies  are  needed 
in  an  1 1  sample  truncation.  A  brief  study  of  the  effect  of  truncation  length  was  conducted  prior  to 
designing  the  hardware  described  in  the  next  section  of  this  report. 

In  principle,  only  the  filter  shown  in  Fig.  3  is  needed.  A  form  of  the  Nyquist  theorem  states 
that  complex  samples  taken  at  a  rate  W  is  sufficient;  hence  we  need  only  values  of  /  and  Q  for  even 
values  of  n.  with  the  value  at  stage  six  being  ±/  and  the  output  value  of  the  FIR  filter  being  ±Q. 

In  Fig.  4,  which  represents  the  hardware  actually  used,  switching  provides  a  pair  of  /  and  Q 
values  for  each  sample.  This  oversampling  of  /  and  Q  economically  provides  additional  redundancy 
for  greater  phase  accuracy. 

CLOCK 

(CYCLE) 


-r(tm  +  1/4  f0  )  SWITCH 


Fig  4  -  DCD  hardware  diagram 

HARDWARE  DESIGN  AND  FABRICATION 

The  FIR  filter  shown  in  Fig.  3  was  fabricated  in  a  compact  form.  Earlier  work  1 1 .4)  dealt  with 
as  many  as  21  filter  stages.  However,  in  the  model  described  here,  only  11  stages  were  actually 
implemented.  This  was  because  the  outer  weights  were  so  small  that  their  effect  on  accuracy  could 
!>e  ignored.  Choice  of  a  12-bit  A/D  posed  the  question  whether  to  use  12-bit  or  16-bit  multipliers. 
The  1 2  -  bit  selection  created  a  need  to  contract  (divide  by  two)  before  the  multipliers  and  expand 
i  multiply  by  two)  after  the  multipliers  because  of  the  summation  operation  preceding  them. 

Figure  5  shows  the  hardware  implementation  of  the  switch  in  Fig.  4.  Signs  for  the  multiplexed 
/.  (J  data  were  obtained  from  the  formulas  given  in  the  previous  section. 

'I  wo  identical  detector*,  were  fabricated  to  provide  two  channels  of  /  and  Q  data  for  the  SEN 
R  \  D  radar  located  at  NRL's  Chesapeake  Bay  Detachment  (CBD).  The  detectors  were  constructed  of 
appro ornately  80  74 1. S  series  transistor-transistor  logic  (TTE)  devices  and  3  complementary  metal 
lode  semiconductor  (CMOS)  (Analog  Devices  1012)  multipliers  with  TTL  I/O  compatibility.  Each 
was  lubricated  on  one  universal-type  wire  wrap  board  (Mupac  202). 
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Fig  5  —  Inversion  and  multiplexing  circuit 


INSTRUMENTATION 

To  test  a  phase  detector  such  as  the  DCD,  a  precise  phase  shifter  is  desired  that  can  vary  the 
input  phase  while  the  output  values  are  being  read.  In  practice  such  phase  shifters,  which  are  needed 
to  measure  the  predicted  0.05°  error  of  the  DCD,  are  not  readily  available.  We  decided  to  use  a 
low-noise  signal  source  with  the  capability  of  exact  frequency  settings  relative  to  the  sampling  fre¬ 
quency.  This  source  provided  a  signal  where  the  phase  change  from  one  sample  to  any  other  sample 
was  precisely  known.  Phase  noise  of  the  signal  source  will  corrupt  the  predicted  values;  since  only 
the  relative  phase  needs  to  be  accurate  and  the  time  between  samples  is  short,  the  phase  error  caused 
by  the  source  was  very  small. 

Figure  6  is  a  diagram  of  the  test  instrumentation.  An  HP3335A  synthesizer  is  locked  to  a  low 
noise,  external  10  MHz  reference  oscillator;  this  reference  oscillator  also  is  used  to  generate  the  5 
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Fig  6  —  Test  instrumentation 


Table  l  —  Measured  /  and  Q  Data  (X  =  Q,  Y  —  f) 
(Signal  level.  0.28  V  rms)* 


Sample 

Y 

Deg.  Calc. 

Deg.  Measure 

Error 

-  125 

13 

-354 

180.00 

177.90 

-2.10 

-  120 

-24 

-346 

-172.80 

-176.03 

-3.23 

-  1 10 

-123 

-329 

-158.40 

-159.50 

-1.10 

-  100 

-  191 

-290 

-144.00 

-146.63 

-2.63 

-90 

-267 

-229 

-129.60 

-130.62 

-1.02 

-80 

-305 

-160 

-1  15.20 

-117.68 

-2.48 

-70 

-345 

-77 

-100.80 

-102.58 

-1.78 

-60 

-344 

7 

-86.40 

-88.83 

-2.43 

-50 

-339 

94 

-72.00 

-74.50 

-2.50 

-40 

-298 

173 

-57.60 

-59.86 

-2.26 

-30 

-249 

246 

-43.20 

-45.35 

-2.15 

-20 

-181 

299 

-28.80 

-31.19 

-2.39 

-10 

-96 

334 

-14.40 

-16.04 

-1.64 

0 

0 

347 

0.00 

0.00 

0.00 

10 

64 

340 

14.40 

10.66 

-3.74 

20 

170 

309 

28.80 

28.82 

.02 

30 

220 

260 

43.20 

40.24 

-2.96 

40 

295 

193 

57.60 

56.81 

-.79 

50 

321 

1  18 

72.00 

69.82 

-2.18 

60 

351 

31 

86.40 

84.95 

-  1.45 

70 

341 

-51 

100.80 

98.51 

-2.29 

80 

322 

-  134 

1 15.20 

1 12.59 

-2.61 

90 

277 

-215 

129.60 

127.82 

-  1.78 

1(H) 

212 

-278 

144.00 

142.67 

-  1.33 

1  10 

141 

-319 

158.40 

156.15 

-2.25 

120 

53 

-342 

172.80 

171.19 

-  1.61 

125 

14 

-344 

180.00 

177.67 

_ 

-2.33 

'Frequency  difference  from  6.25  MH/  is  20. (XX)  H/ 
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\'H  ■.  f...  quern  v  '-t  die  IK'!'  The  HP*33v\  can  he  set  to  a  in  trequenev  in  the  pr«*>  e-sinc 

■  ■ : IX'!')  within  0  * x 1 1  H/.  The  signal  is  digitized  bv  a  12  hit  Attains  Dev  Res  M<  )l)  12t)> 
v  ’ 1  in.'  sv  digitized  values  am  the  inputs  ot  the  IK'D  /  and  <J  i  Hit  puts  ot  the  IK'D  were 

■  .•  •••:.•,•  h»  .i  ]'chi:om\  12  TO  logic  analv/cr 

I  tie  Ins ic  analv/cr  was  set  to  trigger  when  the  Q  value  was  /cm;  then  selected  samples  on  both 
sides  »l  the  trigger  point  were  read  out  Table  1  is  an  example  of  the  data.  The  sample  numbers  are 
m  units  ot  t!ie  5  MU/  sampling  frequency  X  and  Tare  the  Q  and  /  outputs  respective!)  of  the  DCD 
The  Dee  Calx  column  contains  phase  values  based  on  the  offset  trequenev  and  assuming  /ero  at  the 
trisger  sample.  The  Deg  Measure  column  lists  angles  calculated  from  DCD  outputs.  Note  that  bv 
definition.  the  error  is  alvvavs  zero  at  the  trigger  sample:  however  the  actual  phase  ot  the  signal  will 
not  he  exactlv  zero  because  of  A  D  converter  quantization  and  aperture  uncertainty  So  the  measured 
data  will  have  small  bias  errors  that  are  not  due  to  the  DCD  and  thus  are  not  c  tsidered  in  the 
evaluation  ot  DCD  performance  in  the  following  sections. 

KKSIT.TS 

Table  2  contains  DCD  data  taken  bv  using  two  different  A  D  converters  Since  A  D  Si  vielded 
better  results,  it  was  used  throughout  subsequent  tests. 


Table  2  —  A/D  Converter  Comparison 
(Signal  power  =  12  dBm) 


Data  Set* 

A  i  A  / 

0\ 

— 

dr 

Comment 

15  6  020 

25 

10 

.528 

.525 

.507 

A/D  50/Red 

15  6  02  1 

39 

5 

.398 

.398 

.476 

A/D  51 /Blue 

15  6  02  1 

50 

5 

.585 

_ 

.585 

.557 

A/D  51 /Blue 

*  I  Jala  set  information  signifies  date  band  center  signal  off  set' AD  converter: 
.V  is  the  number  of  samples  between  stored  values 


Figure  7  shows  the  effect  of  center  trequenev  on  phase  accuracy.  Note  that  two  phase  accuracy 
measures  are  plotted  against  ce.aer  frequence  /,,  -  (2 M  -  I)  fj 4.  where  M  is  an  integer  and 
/  -  21C  -  5  MHz  In  Fig.  7,  1.25.  5  75.  6.25.  8.75.  and  11.25  MHz  correspond  to  values  of  M 
from  1  to  5.  respectively. 

RMS  phase  errors  r;,.  .  and  a-.  are  calculated  as  follows.  If  <5,  is  the  /th  phase  measurement 
to.  tan  !  (J  // ).  it  can  be  shown  [4|  that  the  variance  of  ct>,  (sine  wave  separated  in  frequency 
from  /,,  b_v  AO  max  be  expressed  bv 


u|  -  (Ao,  -  27rA/An-/2. 

w  here 


Ao  o  .  ;  <i.  . 


A/  is  the  interval  between  samples,  and 


the  bar  indicates  averaging  over  the  range  of  i. 
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fQ.  CENTER  OF  SAMPLED 
FREQUENCY  BAND  (MHz) 

Fig.  7  —  DCD  phase  error  vs  frequence 


It'  signal  frequency  is  constant  but  unknown,  we  may  compute  as  -  (A0,  -  A<£>,)'/2.  As  the 
data  in  Table  2  show,  very  little  difference  existed  between  o\  and  as  because  our  knowledge  of  A/ 
was  precise;  thus  (a]  +  a:)/2  was  plotted  along  with  on  in  Fig.  7. 

Calculation  of  a^  implies  even  more  precise  knowledge  of  signal  frequency.  Assuming 
D|  =  o’|.  we  again  assume  frequency  is  constant  so  that  the  actual  phase  of  the  signal  is 
o',  -  o i  -r-  2 7r( i  -  1 )  A/A/.  In  this  case.  o\  =  (O,  -  0,'r.  Clearly,  if  A/  is  not  accurately 
known,  o',  will  be  in  error  for  large  i.  Note  that  some  difference  exists  between  values  of 
(U|  *  a s ) / 2  and  a\  for  the  higher  frequencies. 

All  of  the  points  plotted  in  Fig.  7  were  plotted  with  A /  =  20  kHz.  Tabic  3  lists  data  sets 
intended  to  demonstrate  the  effect  of  greater  separation  of  signal  frequency  from  band  center.  Since 
sampling  rate  /’  =  5  MH/.  the  Nyquist  band  W  =  2.5  MHz.  so  the  maximum  value  A /  ~  1.25 

MH/  Note  that  the  largest  value  of  A/ in  Table  3  is  0.539  MHz  corresponding  to  one  edge  of  a  cen¬ 
tered  signal  hand  of  w  idth  B  =  1 .078  MHz  Hence,  this  signal  is  oversampled  by  -  1  6 T 

Calculations  were  done  to  predict  errors  that  an  ideal  DCD  would  produce  as  a  function  of  input 

frequency.  The  input  signal  was  sampled,  and  the  samples  were  truncated  to  12-bit  values.  Integer 
arithmetic  was  used  to  simulate  the  calculations  performed  by  the  DCD.  Figure  8  shows  the  mis 
error  as  a  function  of  the  input  frequency;  each  point  is  calculated  from  150  consecutive  outputs  of 
the  DCD  The  center  frequency  is  6.25  MHz.  and  the  sampling  rate  is  5.0  MHz.  Near  center  fre¬ 
quency.  the  predicted  errors  are  -0  05".  errors  remain  below  0.1°  for  an  input  frequency  within 
0  54  MHz  of  the  center,  then  increase  rapidly  for  frequency  offsets  that  are  greater. 

A  new  algorithm  is  being  defined  by  another  investigator  |1()|  that  is  expected  to  further  reduce 
digital  coherent  detector  phase  errors.  At  the  expense  of  another  FIR  filter,  errors  are  expected  to  be 
reduced  by  at  least  an  order  of  magnitude  below  those  plotted  in  Fig.  8. 
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Table  3  —  Accuracy  vs  Frequency  Offset  t A f) 
(Signal  power  =  12  dBm,  A t  =  5,  A/D  51  Blue) 


Data  Set* 

N 

A/(MHz) 

0\ 

°1 

03 

23/6/02/1 

12 

.02 

.209 

.209 

.310 

23/6/15/1 

32 

.150 

.353 

.353 

.308 

23/6/53/1 

12 

.539 

.310 

_ 

.309 

.286 

*Data  set  information  signifies  date/band  center' signal  offset/A/D 
convener.  A /  is  the  number  of  samples  between  stored  values. 


(center  frequency.  6.25  MHz;  sampling  rate.  5.0  MHz) 


SUMMARY  AND  CONCLUSIONS 

We  have  described  hardware  capable  of  converting  radar  and  communications  signals  sampled 
directly  from  RF  or  IF  portions  of  a  receiver.  By  sampling  bandpass  signals  directly,  the  need  for 
dividing  signals  into  parallel  channels  and  converting  to  baseband  is  avoided.  As  a  consequence, 
errors  caused  by  mismatch  between  baseband  channels  are  avoided,  particularly  when  large  received 
power  variations,  typical  of  radar  signals,  are  involved.  Measurements  of  phase  provided  by  this 
hardware  demonstrate  its  intrinsic  accuracy. 

Phase  is  basically  a  measure  of  time  of  arrival  of  the  sine  wave  signal,  relative  to  the  sampling 
pulse  train  (strobe).  If  band  center  is  10  MHz.  <7=1°  corresponds  to  measuring  time  to  an  accuracy 
better  than  I '3  ns.  Note  that  measured  a  <  0.3°  at  /„  ~  3.75  MHz.  which  corresponds  to  time 
accuracy  better  than  230  ps. 
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Time  measurement  accuracies  that  significantly  exceed  the  reciprocal  of  the  signal  bandwidth 
imply  good  signal-to-noise  ratio  (SNR).  Swerling's  estimate  (|9),  pp  4-7,  Eq.  (33)) 
l/C  2  0  SNR)  is  the  expected  rms  time  measurement  error,  and  W  is  the  bandwidth.  The  mea¬ 

sured  value  of  a.  =  UK)  ps  implies  that  SNR  >  34  dB  (IT  =  2.5  MHz).  This  is  well  within  the 
SNR  implied  by  -+•  12  dBm  signal  level  and  a  1  mV  least-significant  bit  (59  dB). 

There  may  be  an  error  minimum  at  3.75  MHz  (Fig.  7).  Phase  noise  spectra  of  most  sine  wave 
sources  diminish  as  frequency  increases.  However,  the  error  caused  by  the  A/D  converter  inceases 
with  increased  frequency  as  the  signal  period  approaches  aperture  uncertainty  time.  The  tentative 
conclusion  is  that  below  3.75  MHz  errors  are  due  mostly  to  phase  noise,  with  AD  aperture  uncer¬ 
tainty  causing  errors  at  the  higher  frequencies. 

In  conclusion,  we  believe  that  even  better  performance  may  be  expected  as  better  direct- 
"ampling  A  D  converters  become  available.  Note  that  calculated  errors  assuming  an  ideal  12-hit  A  D 
converter  (Fig.  X)  are  -0.05°.  compared  with  values  shown  in  Fig.  7.  Further  research  toward 
better  direct-sampling  A  D  converters  is  recommended. 
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